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ABSTRACT 

N') The  distance  travelled  as  a function  of  velocity  for  falling 
spherical  drops  starting  from  rest  was  calculated.  This  was  done  for 
1 to  5 millimetre  spheres  to  determine  the  travel  distance  required  for 
designing  experiments  where  velocities  approaching  terminal  are  re- 
quired. The  results  showed  that  80  percent  of  terminal  velocity  Is 
reached  after  spheres  of  unit  density  and  diameter  of  5 millimetres 
have  fallen  from  rest  through  a distance  of  7.5  metres.  To  achieve  85 
percent  of  terminal  velocity,  an  Increase  of  only  5 percent,  9.5  metres 
Is  required.  Also  calculated  was  the  velocity  of  particles  blown  up- 
ward from  rest  In  a vertical  alrstream.  This  was  done  to  provide  In- 
formation for  comparison  of  experimental  to  theoretical  results  on  the 
effect  of  dust  on  aerodynamic  drag.  ^ 
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simplifying  constants 
drag  coefficient  of  a sphere  In  air 
drag  force  due  to  air  resistance 
spherical  drop  diameter 
acceleration  due  to  gravity 
mass  of  sphere 

effective  mass  of  sphere  In  air  due  to  buoyancy  effects 

Reynolds  number  of  a sphere 

time  after  sphere  starts  to  move 

free  stream  air  velocity 

velocity  of  sphere 

upward  terminal  velocity  of  a sphere  In  a vertical  air  stream 

slip  velocity  of  a sphere  In  a vertical  air  stream 

terminal  velocity  of  a spherical  particle 

displacement  of  sphere  from  rest 

surface  tension  of  the  liquid 

dynamic  viscosity  of  air 
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1.  INTRODUCTION 


The  Immediate  purpose  of  this  report  Is  to  provide  Information 
on  the  velocity  as  a function  of  distance  fallen  from  rest  for  drops  of 
various  sizes  of  Interest,  and  to  describe  the  method  by  which  this  was 
achieved.  The  second  purpose  Is  to  describe  the  determination  of  the 
upward  velocity  of  spherical  dust  particles  in  a vertical  stream  of  air. 
This  Information  was  necessary  for  comparing  experimental  to  theoretical 
results  In  a study  of  Impaction  forces  of  dust  on  spheres  and  cylinders 
(Mellsen), 
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The  mathematical  theory  and  the  solution  of  the  equations 
of  motion,  for  which  computer  programs  are  provided,  are  described 
herein. 

2.  THE  PROBLEM 

It  has  been  shown  experimentally  that  drops  of  various  liquids 
falling  through  air  deform  to  a negligible  extent  provided  that  the 
value  of  pRd2g/Y  is  below  0.4.  For  ordinary  liquids  this  limiting  value 
corresponds  to  drops  which  are  1 - 1.5  mm  In  diameter,  and  the  terminal 
velocities  reached,  under  any  atmospheric  conditions,  by  drops  smaller 
than  this  can  be  calculated  precisely  from  knowledge  of  the  fluid  re- 
sistance offered  to  the  motion  of  spheres  (Batchelor). 

The  work  reported  herein  was  undertaken  to  determine  the  distance 
travelled  to  terminal  velocity  by  droplets  1 - 5 mm  in  diameter  under  any 
atmospheric  conditions.  Although  this  range  of  drop  sizes  falls  outside 
the  limiting  value  described  above,  the  assumption  of  spherical  drops  was 
used.  This  provides  sufficient  accuracy  for  the  purpose  of  designing  ex- 
perimental apparatus.  The  exact  discrepancies  were  not  obtained  but  they 
are  certainly  less  than  the  discrepancies  In  the  terminal  velocities 
themselves,  which  were  obtained  by  comparing  the  calculated  terminal  velo- 
cities of  spheres  to  those  of  water  drops  determined  by  Best's  empirical 
equation  (Batchelor).  The  reason  that  the  discrepancy  in  distance 
travelled  to  terminal  velocity  is  less  than  the  discrepancy  In  the  termi- 
nal velocity  is  that  the  drag  force  gradually  Increases  as  the  drop  travels 
from  rest  to  terminal  velocity,  and,  correspondingly,  the  distortion  of  the 
drop  Increases  to  maximum  deviation  from  spherical  shape  at  terminal  velocity. 

3 . THE  MATHEMATICAL  SOLUTION 

a)  Terminal  Velocity  of  Spherical  Drops 

Equating  the  drag  acting  on  a rigid  sphere  of  diameter  d and 
density  p , falling  with  steady  velocity  vt  to  the  effective  weight  of 
the  sphere  In  air  of  density  p_,  we  have: 

a 
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CD7rd2pavt2  = 7r(p£  " pa)d3g 
8 6 


(Eq.  1) 


where  CD,  the  drag  coefficient,  Is  a function  of  the  Reynolds  number 
Re  = vd/v  which  characterizes  the  motion  of  the  sphere  (Batchelor). 
Hence: 


4(( 


- Pa)gvRe2 


(Eq.  2) 


3paCDRe 


Eq.  2 was  solved  by  an  iterative  technique.  This  was  necessary  because 
the  drag  coefficient  is  dependent  upon  the  Reynolds  number  In  a complex 
way  so  that  the  equation  could  not  be  solved  explicitly  for  vt.  The  de- 
pendency of  drag  coefficient  on  Reynolds  number  has  been  determined  by 
a set  of  definitive  equations  for  CDRe  (Davies)  based  on  the  results  of 
many  experiments.  These  were  solved  for  each  Reynolds  number  as  required 
by  means  of  Newton's  method.  The  initial  estimate  of  CQRe  chosen  for  use 
In  Eq.  2 was  the  Stokes  flow  value  of  24.  This  was  substituted  Into 
Eq.  2 which  was  then  solved  for  v^.  This  value  of  vt  was  then  used  to 
calculate  CpRe  which  was  in  turn  substituted  into  Eq.  2 and  the  process 
repeated.  This  Iterative  procedure  was  tested  and  found  to  converge 
rapidly. 


The  calculations  for  the  terminal  velocity  were  done  by  means 
of  a subroutine  in  the  computer  program  (Subroutine  SBM  36,  Appendix  A). 
The  values  of  CQRe  needed  in  these  calculations  were  obtained  from  a 
function  subprogram  (Function  CDRE  (RE),  Appendix  A). 

b ) Distance  Travelled  from  Rest  in  Still  Air  as  a Function  of  Velocity 

The  equation  of  motion  of  a spherical  particle  falling  in  air 
Is  given  by: 

meg  - D = ^ (Eq.  3) 

WhSre  . ud3(pg  " pa}  (Eq.  4) 

ffle  = 6 
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ird  3p^ 

m * 

(Eq. 

5) 

6 

and 

CDud2p  v2 

D = u a 

(Eq. 

6) 

8 

To  obtain  the  quantity  C^Re  for  use  In  later  calculations, 
the  following  form  for  the  drag  coefficient  was  used: 


c = OWL  (Eq.  7) 

vdp. 

Then,  using  Eq.  7,  substitution  of  Eqs.  4,  5 and  6 Into  Eq.  3 
gives  the  following  equation  for  time  from  rest  as  a function  of  velocity: 


(Eq.  8) 


where 

a ■ g (l  - M 

(Eq.  9) 

l P*j 

and 

b = Jh— 

(Eq.  10) 

SdZ 


Recognition  that  dx  = vdt  gives  the  following  equation  for  dis- 
placement as  a function  of  velocity: 


(Eq.  11) 


Eqs.  8 and  11  were  solved  by  means  of  a computer  program  (Appendix 
A).  The  terminal  velocity  was  first  calculated  by  a subroutine  (SBM  36, 
Appendix  A)  and  the  velocity  increased  by  one  percent  Increments  up  to 
ninety  eight  percent  of  terminal  velocity,  because  terminal  velocity  Is 
asymptotic.  The  quantity  C^Re  was  found  In  each  step  by  the  function  sub- 
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program  (CDRE  (RE),  Appendix  A)  and  substituted  In  Eqs.  8 and  11  with  the 
value  of  v to  solve  for  t and  x at  each  succeeding  velocity. 

c)  Distance  Travelled  Upward  From  Rest  In  a Vertical  Air  Stream 

The  solution  to  this  problem  Is  similar  to  that  of  the  problem 
described  in  the  previous  subsection  with  the  following  exceptions.  The 
upward  terminal  velocity  of  the  spherical  particle  is: 

v*  B u - vt  (Eq.  12) 

where  vt  Is  terminal  velocity  calculated  in  air  at  rest  as  discussed 
earlier.  The  slip  velocity  of  the  particle,  defined  as  the  relative  velo- 
city between  the  air  stream  and  sphere,  Is: 


u - v 


Eq.  6 becomes: 


Vd2pavs 


(Eq.  13) 


(Eq.  14) 


Then  Eqs.  8 and  11  become: 


bvsCDRe 


a - bvsC0Re 


(Eq.  15) 


(Eq.  16) 


The  latter  two  equations  were  solved  by  a computer  program 
(Appendix  B)  which  was  similar  to  the  one  previously  described  (Appen- 
dix A).  The  velocity  was  Incremented  by  one  percent  of  the  upward 
terminal  velocity  Instead  of  one  percent  of  the  terminal  velocity  In 
still  air.  The  final  increment  was  then  at  ninety  percent  of  the  upward 
terminal  velocity,  the  asymptotic  velocity  in  this  particular  case. 
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4.  RESULTS 

The  elapsed  times,  velocities,  distances,  and  fractions  of 
terminal  velocities  were  calculated  for  spheres  of  unit  density  falling 
from  rest  in  still  air.  These  are  shown  for  diameters  of  1 to  5 milli- 
mltres  in  1/2  millimetre  increments  in  the  tables  Immediately  following 
the  listing  of  the  computer  program  by  means  of  which  they  were  cal- 
culated (Appendix  A).  Also  shown  are  the  terminal  velocities  and,  for 
comparison,  the  terminal  velocities  of  water  drops  from  Best's  equation 
(Batchelor).  The  results  are  illustrated  in  the  form  of  velocity  versus 
distance  curves  for  1,2,  3,  4 and  5 millimetre  spheres  (Fig.  1),  and  in 
terminal  velocity  versus  diameter  curves  for  rigid  spheres  compared  to 
water  drops  (Fig.  2). 

The  elapsed  times,  velocities,  distances  travelled,  and  fractions 
of  upward  terminal  velocities  were  calculated  for  spheres  of  specific 
gravity  2.6,  corresponding  to  the  density  of  glass  beads  used  in  experi- 
ments (Mellsen)  on  spheres  starting  from  rest  in  a vertical  air  stream. 

An  example  of  the  results  of  these  calculations  is  shown  in  the  table 
following  the  listing  of  the  computer  program  by  means  of  which  they  were 
obtained  (Appendix  B).  The  results  are  illustrated  in  the  form  of  curves 
of  upward  sphere  velocity  versus  upward  air  velocity  for  three  particle 
diameters  (Fig.  3).  These  diameters  and  the  distance  travelled  were  used 
in  previous  experiments  to  measure  drag  on  spheres  and  cylinders  in  dusty 
air.  The  results  shown  (Fig.  3)  were  used  for  comparison  of  theory  to 
the  results  of  these  experiments  and  reported  elsewhere  (Mellsen).  The 
program  listed  (Appendix  B)  was  also  used  to  provide  direct  numerical 
values  for  that  work. 

5.  DISCUSSION 

As  can  be  seen  (Fig.  1)  better  than  80  percent  of  terminal 
velocity  was  reached  by  the  time  spheres  of  unit  density  and  maximum  dia- 
meter of  5 millimetres  had  fallen  from  rest  through  a distance  of  7.5  metres. 
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To  achieve  85  percent  of  terminal  velocity,  an  Increase  of  only  5 percent, 
a distance  of  9.5  metres  was  required.  This  Indicated  that,  If  a drop 
tube  Is  to  be  used  in  experiments,  a great  deal  of  increase  in  length  is 
required  to  gain  little  increase  in  velocity  when  approaching  terminal 
velocity.  However,  a useful  consequence  of  this  asymptotic  behavior  is 
that  the  percentage  of  terminal  velocity  acqulr  : by  liquid  drops  Is 
closely  approximated  by  calculated  percentages  r spherical  particles. 

The  difference  between  terminal  velocity  for  spheres  of  unit  specific 
gravity  and  water  drops  Increases  with  drop  size  (Fig.  2).  However,  the 
difference  In  percentage  of  terminal  velocity  reached  for  drop  sizes  of  1 
to  5 millimetres  is  still  negligible  from  the  point  of  view  of  Impaction 
test  accuracy. 

The  curves  of  sphere  velocity  versus  upward  air  velocity  for  a 
travel  distance  of  95  centimetres  (Fig.  3)  showed  that  terminal  velocities 
were  not  reached  for  air  velocities  above  480  and  550  centimetres  per 
second  for  particle  diameters  of  0.0470  and  0.0155  centimetres,  respectively. 
Consequently,  terminal  velocities  were  not  reached  in  experiments  to  which 
other  calculations  were  compared  (Mellsen).  The  applications  and  con- 
sequences of  these  results  are  described  In  the  same  reference. 

6.  CONCLUSIONS 

Results  of  velocity  versus  distance  for  drops  falling  in  still 
air  and  spherical  dust  particles  blown  upward  in  an  air  stream  were 
obtained.  These  results  are  useful  in  the  design  or  analysis  of  experi- 
ments where  optimum  travel  distances  are  required  or  where  travel  dis- 
tances are  confined  to  available  working  space. 
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Figure  1:  Velocity  vs  Distance  for  Spheres  of  Unit  Density  Falling  froa  Rest  in  Still  Air  at  15°C 

and  76  c*  Hg  Pressure 
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THE  TERMINAL  VELOCITY  OF  A LIQUID  DROP  BY  BEST'S  EQUATION  IS  530*606  CK/SEC 
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TERMINAL  VELOCITY  OF  A LIQUID  DROP  BY  BEST"S  EQUATION  15  791.973  CM/SEC 
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TERMINAL  VELOCITY  OF  A LIQUID  DROP  BY  BEST *S  EQUATION  IS  691.9**  CM/SEC 
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reached  after  spheres  of  unit  density  and  diameter  of  5 millimetres 
have  fallen  from  rest  through  a distance  of  7.5  metres.  To  achieve  85 
percent  of  terminal  velocity,  an  Increase  of  only  5 percent,  S.5  metres 
Is  required.  Also  calculated  was  the  velocity  of  particles  blown  up- 
ward from  rest  In  a vertical  alrstream.  This  was  done  to  provide  In- 
formation for  comparison  of  experimental  to  theoretical  results  on  the 
effect  of  dust  on  aerodynamic  drag. 
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